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A  FAMILY  OF  CURVES  FOR  THE  ROUGH  SURFACE 
REFLECTION  COEFFICIENT 

INTRODUCTION 

Recently,  Miller,  Brown,  and  Vegh  [1]  gave  a  short  derivation  that  showed  the  rough  surface 
reflection  coefficient  R,  or  roughness  factor,  was  given  by 

R  =  I  E/E0  r0  i  -  exp  (  — 2(2irg)2  )/0(  2(27rg)2 )  (1) 

where 

g  =  (it  sin  0)/A 

Here  E  is  the  average  electric  field,  E0  the  field  due  to  the  direct  wave,  T0  the  smooth  sea  reflection 
coefficient,  o-  the  standard  deviation  of  the  sea  surface  elevation,  0  the  grazing  angle  for  specular 
reflection,  A  the  electromagnetic  wavelength,  Iq(z)  the  modified  Bessel  function  J0(iz),  and  E/E0r0 
the  normalized  coherent  reflected  field;  g  is  a  measure  of  the  effective  surface  roughness  or  simply  sur¬ 
face  roughness  [2,  p.  10]. 

Equation  (1)  was  first  obtained  by  Brown  and  Miller  13]  in  1974.  This  result  agrees  well  with  the 
experimental  results  obtained  by  C.I.  Beard  [4]  in  the  range  0  <  g  <  0.3  rad.  An  earlier  theoretical 
result  of  W.S.  Ament  [5]  had  agreed  with  Beard’s  experimental  curve  only  in  the  range  0  <  g  <  0.1 
rad.  Ament’s  result  is  simply  Eq.  (1)  without  the  70  factor  (Fig.  1). 


Fig.  1— Comparison  of  theoretical  and  experimental  results 
for  the  reflection  coefficient 


Equation  (1)  was  derived  in  Ref.  1  by  taking  the  Fourier  transform 

E/E0  r0  -  J  b  exp  |  4  ni  sin  0  j  Diy)dy  (2) 
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of  the  probability  density  function  D(y)  for  sea  wave  elevation  y  given  by  Eq.  (3).  Equation  (2)  was 
obtained  by  Ament  [51  and  may  be  understood  by  observing  that  D(y)dy  is  the  frequency  of 
occurrence  of  a  plane  wave  that  is  reflected  from  the  sea  surface  with  a  grazing  angle  ip  and  at  a  point 
between  y  and  y  +  dy.  The  expression  (4  ny  sin  0)/A.  represents  the  phase  of  the  electromagnetic 
wave  after  reflection. 

The  statistical  model  for  sea  wave  elevation  y  used  in  obtaining  Eq.  (1)  is 

y  *  H  sin  x  (3) 

A 

where  sea  wave  crest  height  H  is  distributed  normally  and  x  is  distributed  uniformly  in  the  interval 
I — A/4,  A/4],  A,  denoting  water  wavelength.  This  model  implies  that  the  sea  surface  is  divided  into  a 
large  number  of  domains  and  that  within  each  domain  the  elevation  variation  is  a  single  sinusoid  with 
random  Guassian  amplitude.  Other  statistical  models  of  the  ocean  surface  may  be  found  in  Refs.  6  and 
7. 


In  this  report,  we  use  a  theoretical  density  for  H  derived  by  Rice  [8]  and  by  Cartwright  and 
Longuet-Higgins  [9],  together  with  the  model  of  the  sea  surface  elevation  given  by  Eq.  (3)  to  obtain  a 
family  of  curves  for  the  rough  surface  reflection  coefficient.  We  obtain  an  envelope  wherin  these 
curves  lie.  Additionally,  we  obtain  the  distribution  function  for  sea  wave  elevation  under  the  stated 
assumptions. 


PRELIMINARY  RESULTS 

Rice  and  Cartwright  and  Longuet-Higgins  have  derived  an  expression  for  the  probability  density 
of  H  which  may  be  expressed  in  the  form 


K  (//,«)  = 


<T  H'j2n 


exp 


- H 2 


2  «2o-h 


-H1 

2<t2h 


1  +  erf  ~ 

2  <th 


y/2  H  yll  ~  cJ 


(4) 


Here  erf  (z)  is  the  error  function  defined  by 

erf  (z)  =  JQ  exp  (  -i2  )  dt 

o’//  is  the  standard  deviation  of  F,  and  0  <  e  <  1.  The  parameter  e,  known  as  the  spectral  width 
parameter,  is  described  fully  in  Refs.  9  and  10,  p.  515.  It  may  easily  be  seen  that 
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In  the  first  instance  we  have  a  Gaussian  density  as  a  limiting  case;  in  the  second  instance  we  have  a 
Rayleigh  density  (Fig.  2). 


Fig.  2— Density  function  for  the  sea  wave  crest  height  H  for  different 
values  of  the  spectral  width  parameter  e 


It  is  shown  in  Appendix  A,  with  the  model  of  sea  wave  elevation  given  by  Eq.  (3)  and  the  density 
function  for  wave  crest  height  H  given  by  Eq.  (4),  that  the  standard  deviation  of  sea  wave  elevation  <r 
is  given  by 

<r  —  crHf{  V5  t}  ) 

where  •»)  is  defined  by 


V  =  1  +  y  (  I  -  e2  ) 


Then  it  is  shown  in  Appendix  B  that  the  probability  density  function  for  y  is  given  by 


.  Vi  —  c 

+ - exp 


^  Mo 


exp  (-s2)  erf  -  —  -  \s2  + 


v)  | 


where  A"0  is  a  Bessel  function  of  imaginary  argument  defined  by  Eq.  (D5)  in  Appendix  D.  Graphs  for 
the  cumulative  distribution  function  for  sea  wave  elevation 


Piy,t)  D(t,e)  dt 


are  given  in  Fig.  3.  An  alternate  form  for  Eq.  (7)  is  developed  in  Appendix  D.  When  substituted  into 
Eq.  (8),  a  formula  for  P(y,c)  is  obtained  that  is  suitable  for  numerical  quadrature.  This  formula  for 
P(y,t)  is  also  given  in  Appendix  D.  Yet  another  form  for  Eq.  (7)  is  given  in  Appendix  E,  together 
with  a  formula  for  the  special  case  D(y,yfl!D  which  is  of  interest  and  will  be  discussed. 
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Fig.  3— Cumulative  distribution  function  for  sea  wave  elevation  y 
for  different  values  of  the  spectral  width  parameter  c 

THE  REFLECTION  COEFFICIENT 

Substituting  the  expression  for  D(y,t)  given  by  Eq.  (7)  into  Eq.  (2),  we  obtain  the  following 
expression  for  the  rough  surface  reflection  coefficient  R : 

R  (#,«)  =  «2  exp  (  -2€2tj2  (2 ng)1  )/0(  2€2t>2  (2 ng)2  )  +  Vl  -  e2  exp  (  -4r»2  ( 2irg )2  ) 

-  |  (1  -  e2)  $1  y,  1;  2;  e2, -4e2i72  (27Tg)2  (9) 

Here  4»|(a,  >3;  y;  x,  y)  is  a  degenerate  hypergeometric  function  in  two  variables  x  and  y,  defined  for 
|x|  <  1,  |y|  <  oo  by 

«&,(«,  /3;  y;  x,  y)  =  £  (a,n_+m)(frn)  (I0) 

„.„.o  (y'n  +  m> n  m- 

where  the  Pochhammer  symbol  (a,n)  is  defined  by  a  ratio  of  gamma  functions: 

( a,n )  =  T(a  +  n)/r(a)  -  a(a  +  1)  ...  (a  +  n  -  1) 

This  result  for  the  reflection  coefficient  is  obtained  in  Appendix  C. 

Graphs  for  the  reflection  coefficient  R  (g,t)  as  a  function  of  surface  roughness  g  are  given  in  Fig. 
4.  The  envelope  of  these  curves  are  effectively  which  forms  the  upper  envelope,  and  /?(g,0) 

which  forms  the  lower  envelope.  The  upper  envelope  corresponds  to  the  Guassian  density  given  by  Eq. 
(5)  and  for  which  the  reflection  coefficient  is  given  by  Eq.  (1) 

/?(*,!)  -  exp  (  -King)2  )/„(  2(2 ng)1  ) 
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Fig.  4— Theoretical  curves  for  the  reflection  coefficient  for  different  values 
of  the  spectral  width  parameter  t 


The  lower  envelope  corresponds  to  the  Rayleigh  density  given  by  Eq.  (6),  and  for  which  the  reflection 
coefficient  is  given  by 

*(g,0)  -  exp  J  2  (2ng)2 

■  -4.: 

Both  these  limiting  envelopes  may  be  incorporated  into  the  general  equation,  provided  we  define  the 
third  term  in  Eq.  (9)  to  be  zero  when  «  -  1. 


SOME  OBSERVATIONS 

Note  that  for  surface  roughness  0  <  g  <  0.1  rad.  Rig*)  for  all  «  and  Ament’s  result 
Rig)  -  exp  l-2(2wg)2)  are  similar  and  agree  with  Beard’s  experimental  curve.  This  is  not  surprising 
since  all  theoretical  results  must  be  consistent  with  the  Rayleigh  criterion  12]  On  the  interval 
0  1  <  g  <  0.3  rad,  numerical  calculations  reveal  that  the  experimental  curve  of  Beard  corresponds 
almost  exactly  to  *(*,0.7).  Noting  that  0.7  =  75/2,  a  closed  form  for  Diy, 72/2)  ts  derived  in 
Appendix  E,  since  this  case  is  of  special  interest  (Fig.  5).  Other  data  points  in  the  domain 
O.i  <  g  ^  0.3  rad  given  by  Beard  [4,  Fig.  2]  lie  within  the  two  envelopes  *  (*,1)  and  *(*,0)  (Fig.  W. 

CONOLU^^,^  ^  (.;c/  ^ 

Statistical  models  for  the  sea  surface  hav*been^ut  fefw^amfcused  to  derive  a  family  of/curves 
for  the  rough  surface  reflection  coefficient.  One  of  these  curves  fits  the  experimental  data  curve. 
Other  experimental  data  points  lie  within  the  envelope  of  the  derived  curves.  There  is,  therefore,  a 
theoretical  basis  for  the  experimental  data  when  the  surface  roughness  exceeds  0.1  rad.  /f  4  s  / 


Fig.  5— Comparison  of  Beard’s  experimental  curve  for  the  reflection 
coefficient  and  the  theoretical  curve  for  e  —  1/V2 
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Appendix  A 

THE  STANDARD  DEVIATION  OF  SEA  WAVE  ELEVATION 


Consider  the  random  variable  y  =  H  sin  6,  where  H  is  a  random  variable  with  density  K(H, c) 
given  by  Eq.  (4),  and  0  is  a  random  variable,  independent  of  //,  distributed  uniformly  on  the  interval 
|0|<7r/2.  Let  (J(0)  be  the  density  of  0,  so  that 


U(9)  =  it  1  \o\  ^  ?r/2 

=  0  \o\  >  -nil 

The  expected  value  or  mean  of  a  random  variable  x  is  defined  by 


Hx  =  E[x  1  =  f  x  /Gc)  dx 

where  /  is  the  density  function  for  x  and  the  variance  is  defined  by 

crl  =  £[  Ox  ~E[x l)2  ]  =£[  (x  ~nx)2  ]  =£[x2]  -fi2 


If  x  and  y  are  independent,  then 


£Uy]  =  £[x]  £[y] 


In  our  case  since  H  and  sin  0  are  also  independent. 


C 

tiy  =  E[H  sin  0]  =  E[H]  I  ,  sin  0  U(0)  dO  =  0 

u  —  7T/2 

a-2  =  E[H2  sin2  0]  =  E[H2]  f  '  sin2  0  U(0)  d0  =  E[H2] 

—nil  7 


Hence 


2<r2  -  <r2H  +  n}, 


-  /_”//£(//,€>  dH 


p  44  « 


2V2  o-w  Vl  -  €2  f  x2  exp  (-x2)  -  Vw/2  crH  Vl  -  e2 


Appendix  B 

THE  PROBABILITY  DENSITY  FUNCTION  FOR  y 


Consider  the  random  variable  y  -  H  sin  0,  where  H  is  a  random  variable  with  density 
K{H,e),  |//|  <  «>  given  by  Eq.  (4),  and  0  is  a  random  variable,  independent  of  H,  with  density  U(9) 
given  by  Eq.  (Al).  Let  Z)(y,«)  be  the  density  function  for  y.  Since  0  and  H  are  independent,  the 
joint  density  g(9,H ,€,)  of  9  and  H  is  given  by 


g(9,H,f)  -  K(H,t)  U(9) 
0 


|0|  <  tt/2,  i//|  <  oo 

elsewhere 


Set  v  —  H  so  that  9  —  sin  (y/v)  and  let  /(y,v,e)  be  the  joint  density  of  y  and  v.  Then 


/(y,v,«)  -  |/|  g(sin  1  (y/v),  v,e) 


where  the  Jacobian  determinant 


A  straightforward  computation  yields 


I.,  a 9  a h  a 9  &h 

dy  dv  dv  dy 


Ul  =  (v2  —  y2)~l/2 


and 


/(y,//,€)  -  n'1  (H2  -  y2)~1'2  K(H,t) 
0 


|y|  <  \H\  <  00 
elsewhere 


Now,  integrate  /(y,//,c)  over  H  to  obtain 


Dty  f)  =  jl  rw  KjM, dH  +  -i  r  K(H*} 

v  7T  J—  yfWZp  7T  y/W^ 
Substituting  K(H,t)  into  this  equation  and  using  Eq.  (A2)  we  have 


dH 


D(y'e)~  zikz;  f\ 


exp 


-H2 

dH  pyn 

-H2 

erf 

4c27)2<r2 

Vl  -  e2  f 00 

V°-2 

f  2-qa-e 

ir3/27j(r  y/H2  —  y2 


dH 


2tT7)2<T2 


y/H2  -  y2 


where  c r  and  rj  are  defined  by  Eq.  (A3).  A  change  of  variables,  namely,  H  —  [y|  y/x  in  the  former 
integral  and  H2  -  y2  —  s2  in  the  latter  integral,  gives 


10 


Appendix  C 

THE  REFLECTION  COEFFICIENT  *(*,«) 


Substituting  the  expression  for  D(y,e )  given  by  Eq.  (7)  into  Eq.  (2),  we  obtain 

+ x;  x'-p  <-*’> 


cos  y  sin  ifi  erf 


-«*  si  i  y 2 

<  4tj2o-2 


I J/2 


ds  dy 


The  first  integral  is  evaluated  by  making  the  transformation  x  -  y2/ 8e2T»2a2  and  then  using  Ref.  11,  p. 
765,  6.755-9.  See  Ref.  1  for  essentially  the  same  computation.  We  have  then 


9-/l  _  -1  (•  oo  C  °o 

R  -  e2  exp  (  -2«  V  (27rg)2  )/0(  2e  V  (2irg)2  )  +  — — — —  J  J  (  -  )  ds  dy 


nr)cr 


where  g  is  defined  by  g  =  (cr  sin  i|i)/A  and  (•)  is  the  integrand  of  the  double  integral.  Making  the 
transformation  t  -  y/2-q<r  and  then  changing  to  polar  coordinates  gives 


2-v/l -«2  r°° 


f  f  (  )ds  dy-  f  f  cos  (Sirrigt)  exp  (- f2 -  s2)  erf 

•'0  7T  •/0  •'0 


iT7)(r 


y/l  —  e2 


>/r2+s2 


ds  dr 


*  J*o  ''  J*o  rexp(-r2)cos(87TT)gr  cosfllerf  — ^  —  r  dr 


d0 


Since  by  Ref.  11,  p.  402,  3.715-19, 


J  cos  (iirrigr  cos  9)  dO  —  -y  J0(iir^gr) 


we  have 


R(g,t)  -  «2  exp  (  -2e2rf2  (2i rg)2  )/o(  2«2tj2  (27rg)2  ) 


+  2Vl  -  €2  f  r  exp  (-r2)  J0dm)gr)  erf  — - - —  r 

0  € 


(Cl) 
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To  obtain  a  closed  form  for  R ,  we  show  the  following:  for  any  real  numbers  a  and  p  ^  0, 


J*  r  exp  (-r2)  J0(ar)  erf  (fir)  dr 


—  exp  — —  I  fl(l  4-  P2)~i/2  4>i  —  1-2’  — - —  - — — 

2  CXP  4  )  *[  2’1,A  1  +  p2'  4(1  +fi2) 


where  <J>!  is  a  degenerate  hypergeometric  function  of  two  variables  defined  by  Eq.  (10).  Substituting 
Ref.  11,  p.  931,  8.252-4 


erf  (fir)  -  1  -  &  exp  (~fi2r2)  J~ 


00  exp  (— r2t 2)  dt 


2  +  B2 


r  >  0 


in  the  left  side  of  Eq.  (C2)  and  noting  that  Ref.  11,  p.  717,  6.631-4  for  p  >  0  and  any  real  q 


C00  1 

J  r  exp  (-pr2)  J0(qr)  dr  =  —  exp 


we  obtain 


o  r  exp  [  -(1  +  fi2  +  t 2)  r2  ]  J0(ar)  dr  J  +  ^ 


Now  applying  Eq.  (C3)  again  to  the  integral  in  braces  we  obtain 


1  ZbL  _  JL  f  ~  ‘  1  4(1  +  p2  +  t2)  j 

2  CXP  4  n  (l  +  b2  +  t2)  (fi 2  +  t2) 


and  on  making  the  transformation  t  -  Vl  +  /32  Vjc-1  -  1  we  have 


_1_  — a 2 _ £ _ 

2  CXP  4  2n  (1  +  fi2)*2 


f '  x1'2  (1  -  x)"1'2  [l  -  -Hh-I  '  «P 

J°  1  +  /32  j  4(1  +  p2) 


Finally,  using  Ref.  11,  p.  321,  3.385  we  obtain  Eq.  (C2). 


Letting  a  —  iir-qg,  p  =  Vl  -  e2/e  in  Eq.  (C2)  and  substituting  the  result  into  Eq.  (Cl)  gives  Eq. 
(9).  We  may  also  write 


R  (g,e)  -  e2  exp  (  -2e27j2  (2ng)2  )/0(  2e2Tj2  (2ng)2 )  +  Vl  -  e2  exp  (  -4tj2(2 ng)2 ) 


-  it  '«2  (1  -  e2)  J  x,/2  (1  -  x)~,/2  (1  -  e2x)-1  exp  (  -4«2i72(2w^)2x  )  dx  (C4) 


This  expression  for  R  is  used  in  Appendix  D. 


mmmi 


r 


Appendix  D 

THE  CUMULATIVE  DISTRIBUTION  P(y,e) 


Defining  m  =  (4 tt/X)  sin  0,  we  obtain  from  Eq.  (C4) 

R(<o,e )  —  e2  exp  (  — y  e2r)2cr2at2  )Iq(  y  e2i? 2<r2<o2  )  +  Vl  —  «2  exp  (  —  tj2ct2co2 ) 

- —  e2  (1  -  «2)  f  xl/2  (1  —  x)~l/2  (1  —  «2x)-1  exp  (— e2ij2cr2co2x)  dx 

„  Jo 

Since  we  may  write  Eq.  (2) 

R(co,e)  -  f“  e‘*y  D(y,e)  dy 

•/  —  oo 

we  must  have 

D(y,e)  =  ——  f  e~iuy  R(w,e)  dot 
Substituting  the  expression  for  R  (w,«)  into  this  equation  gives 

D(y,e)  -  — A - exp  (-y2/&€2r,2o-2)K0(y2/&e2r)2<r2)  +  —  Vl  -  c2  J  cos  (coy)  exp  (-t >Vco2)  do> 

2  It3  1)(T  ^ 

_  i2il  ~  J*o  cos  (<oy)  |  j*0  x1/2  (1  -  x)_1/2  (1  -  c2x)_1  exp  (-«Va2co2x)  |  du> 

Noting  that  Ref.  11,  p.  496,  3.952-9 

f  exp  (-/32x2)  cos  (ax)  dx  =  exp  (  -a2/ 4/32 ) 

•'O  zp 

we  may  integrate  the  single  and  double  integrals  with  respect  to  co  to  obtain 

Vl  -  € 2  (  -y2  I  e(l  -  «2)  f1  exp  (-y2/4t2n2<r2x) 

2irl/V  eXP  j  4t}2o-2  j  2tt3/V  Jo  0  ~  «2*)  VT=^ 

for  the  latter  two  terms  in  the  equation  for  D(y,c).  On  making  the  transformation  x(r  +  1)  <■  1,  we 
have  another  form  for  D(y,«),  viz. 
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Now,  making  the  transformation  cosh  x  -  2t  +  1  and  using  Ref.  11,  p.  289,  3.222-2,  we  have 


Finally,  combining  equations  (D2),  (D3),  (D4),  (D6),  and  (D7),  we  obtain 


r 

Jo 


dx 


1 


VI  +  cosh  x  V2  V7  (t  +  1) 


f 

Jo 


dt 


it 

V2 


(D7) 


1  .  Vl  -  e2  , 

e2(l  —  e2)  i 

erf 

r“ 

y  cosh  s 

2ei7<r 

ds 

.  +  _  en 

2  2 

2tjo- 

*  J 

0  cosh  s  (cosh2  s  -  e2) 

.  vVl  +  cosh  s  . 

erf  ~V  ns -  ds 

V2e2  r°°  2V26T)cr _ 

+  2 ir  •'O  VI  +  cosh  s 

This  expression  for  P(y,e),  which  is  suited  for  numerical  evaluation,  was  used  to  obtain  the  graphs  in 
Fig.  3. 
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Appendix  E 

A  CLOSED  FORM  FOR  THE  DENSITY  FUNCTION  Diy<j2j2) 


To  give  a  dosed  form  for  the  probability  density  function  D(y,V2/2),  we  first  show  that  for  any 
real  k  and  n 

n{/2  fQ  exp  (—s2)  erf  (k  yJfi2  +  s2)  ds 


lr  r  £m2<1  +  *2)  i  1  —  k2  1 

- t*0"1  k+  YTic2  Jo  exp  [  JTk2  5  ]  Ko(s)  *  (E1) 

To  show  this,  note  that 

-j-  erf  (k-Jn2  +  s2)  —  — -■  exp  (- k2n 2)  exp  (- k2s 2) 

dfi  vir  y/fi2  +  s 2 

and  hence 


erf  (k  yJpf+  s2)  -  exp  ( ~k2s 2)  J*  *  1  exP  ^  kJ  ^  dt  +  erf  (fc|sl) 

Vw  Vs2  +  r2 

Making  the  transformation  /  -  \ti\yfy  gives 

erf  ( k  y/fi2  +  s2)  -  exp  (.-k2s2)  f  eXff-y--*V^  dy  +  erf  (fc|s|) 
Multiplying  this  by  exp  (— s2)  and  integrating,  we  have 

J*0  exp  (-s2)  erf  ( k  yin2  +  s2)  ds  -  -^=-  JQ  exp  (~k2n2y)  |  /Q  ?xp  -  —  ds  J  <(y 

+  J*Q  exp  (-s2)  erf  (its)  ds 


Using  Ref.  11,  p.  649,  6.285-1,  we  have 


J0  exp  (-s2)  erf  (its)  ds  - 


tan ~lk 
y/n 


The  integral  in  braces  is  evaluated  by  making  the  transformation  s2  -  x  and  then  using  Ref.  11,  p.  316, 
3.364-3  to  obtain 


MILLER  AND  VEGH 


J*Q  exp  (— s2)  erf  (k  V/x2  +  s2)  ds 


tan 


rlk  ku.2  r1 

K~+  itk  Jo  exp 


y/i2  (1  -  *2)  y 


j/jl2  (1  +  k 2)  y 


Vir  2Vtt 

Finally,  making  the  transformation  m2(1  +  k2)  y  -  2s,  we  have  Eq.  (El). 

Using  k  -  V 1  -  e2/e,  n  -  y/2r)cr  in  Eq.  (El)  and  substituting  the  result  into  Eq.  (7)  gives 


D(y,e) 


-v2 


2.3/  V  eXP  I 

1  8e2i)2a- 

Vl  —  e2 

- y 2 

3/2  CXP 

7r  '  r/er 

4ij  2<r2 

8e2ij2<r2 


+  -^.y— -  exp  — 4—z  ]  (  cos-1€  +  eVl  —  e2  f; /8*  v  a  exp  [  (2c2  —  l)s  ]  K0(s)ds  1  (E2) 

ir1/zi)(r  4r)  <rz  J  l  0  ) 


It  may  be  shown  by  using  Ref.  12,  p.  87  that 

J*o  Xo(s)  ds  -  y-x[  K0(x)  L_j(x)  +  K\(x)  L0(x)  ] 

where  Af„(x)  are  modified  Bessel  functions  and  L„(x)  are  modified  Struve  functions.  Defining 

a  s  4/V4  +  7T, 


we  obtain 


Z>(y,V2/2) 


exp  ( — y2/a2cr2 ) 
2(2ir)1/2  aa 


(E3) 


which  is  a  closed  form  for  D(y,V 2/2). 

The  functions  f  K0  (/)  dt,  Kv(x),  L„(x)  may  be  computed  using  tables  found  in  Ref.  13  so 
that  either  Eq.  (E2)  or°Eq.  (E3)  may  be  used  to  compute  Z)(y,V2/2).  We  remark  that  it  may  be  shown 
for  0  <  |a  I  <  1,  0  <  x 


Jo  exp  (at)K0(t)  dt 


sgn  a 


exp  (ox) 


cos  (a  x  /)  dt 

(1  +  f2)  Vl  +  a2t2 


■v'U5&r--n 


END 


i  I 
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